Statistical Properties of Random Banded Matrices with Strongly 



Fluctuating Diagonal Elements 



Yan V. Fyodorov 1 '* and Alexander D. Mirlin 2 '* 

1 Fachbereich Physik, Universitdt-Gesamthochschule Essen, Essen 45117, Germany 
2 Institut fur Theorie der Kondensierten Materie, Universitdt Karlsruhe, 76128 Karlsruhe, 

Germany 
(February 1, 2008) 

Abstract 

The random banded matrices (RBM) whose diagonal elements fluctuate much 
stronger than the off-diagonal ones were introduced recently by Shepelyansky 
as a convenient model for coherent propagation of two interacting particles in 
a random potential. We treat the problem analytically by using the mapping 
onto the same supersymmetric nonlinear a— model that appeared earlier in 
consideration of the standard RBM ensemble, but with renormalized param- 
eters. A Lorentzian form of the local density of states and a two-scale spatial 
structure of the eigenfunctions revealed recently by Jacquod and Shepelyan- 
sky are confirmed by direct calculation of the distribution of eigenfunction 
components. 
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The ensemble of Random Banded Matrices (RBM) can be generally described as that 
of large N x N matrices having nonzero elements effectively within some wide band of the 
width b ^> 1 around the main diagonal. Such a structure naturally appears in various 
physical contexts, and serve as a useful model in Quantum Chaos , atomic physics and 
solid state physics || . Due to this fact a lot of efforts were spent in order to study different 
kinds of RBM, both numerically 0|| and analytically f||7]]. In particular, it was found that 
the problem can be mapped onto a supersymmetric Id nonlinear cr-model introduced in || , 
provided all matrix elements within the band are independent and distributed around zero. 
More precisely, the mapping was shown to exist for those matrices whose variance (\H, 



|2\ 



was dependent on the distance \i — j\ from the main diagonal: (\Hij\ 2 ) = b~ l f(\i — j\/b) 
where the function f(r) is of the order of unity when r < 1 and decreases exponentially (or 
faster) at r 3> 1. 

Quite recently, Shepelyansky |J argued that a very interesting problem of two interact- 
ing particles propagating in a quenched random potential can be effectively mapped onto a 
class of RBM whose diagonal elements Ha fluctuate much stronger than off-diagonal ones: 
(\Ha\ 2 ) / (\Hij\ 2 ) oc 6 3> 1. Using this kind of mapping, Shepelyansky predicted a consider- 
able interaction-assistant enhancement of the two-particle localization length as compared 
with the localization length of one particle in the same random potential. This conclusion 
was confirmed later on by Imry [ITJ who employed the Thouless scaling block picture bypass- 



ing the mapping to RBM. Subsequent numerical studies |TlJ also confirmed main qualitative 
result by Shepelyansky, but revealed some deviations from the predicted behavior of the two- 
particle localization length, which were attributed to oversimplified statistical assumptions 
concerning RBM elements in Shepelyansky construction. Nevertheless, it is clear that She- 
pelyansky RBM model (SRBM) catches adequately at least some of the important features 
of the original physical problem and thus deserves more detailed study. 

In a very recent paper [fT2] , Jacquod and Shepelyansky presented their detailed nu- 
merical results on statistical properties of SRBM. They revealed a peculiar structure of 
eigenfunctions consisting of a set of large spikes separated by regions of relatively 
small amplitude. Such a "sparse" spatial arrangement shows up in a difference between 
the localization length I related to the rate of a spatial decay of an eigenfunction enve- 
lope, I = - lim^oo In | 1 I r Q ,(0) 1 i r a ,(n)|, and the length £ defined as the participation ratio, 
£ = (J2 n \^a(n) | 4 ) . For a conventional "dense" eigenfunctions these two lengths are ex- 
pected to be of the same order of magnitude, whereas for SRBM it was found that I £. 
Another interesting feature making SRBM different from earlier studied cases is that in any 
given realisation of the disorder the local density of states (LDOS) defined as 

p(E,n) = J2 I *a(«) \ 2 5(E-E a ) (1) 

a 

was found to follow the simple Lorentzian form with a width T oc W^ 1 independent of the 
parameter b, where Wb 3> 1 determines the scale of fluctuations of the diagonal elements 
H nn . To this end it is appropriate to mention that the Lorentzian form of LDOS was earlier 
found to be typical for RBM with linearly increasing mean value of the diagonal elements: 
(H nn ) = (3n @,|. 



In the present article we show that Shepelyansky RBM model can be again mapped 
onto the standard one-dimensional nonlinear a— model with modified parameters. This fact 
allows us to reproduce analytically most of peculiar features of the SRBM discussed above. 
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We consider the random Hermitian Ullfl matrix = Widij + , where the matrix 



is a standard RBM characterized via the variances: = (H$*H-j ) = \f(\ i — j \ /&), 
normalized in such a way that Y^T=-oo\f( r /fy — 1- This normalization ensures that the 
width of the energy spectrum of the matrix H-j is of order unity in the limit b — ► oo. The 
parameters Wi are assumed to be independently distributed around zero according to the 
probability density V(W) = ^r b h(W/W b ) where h(r ~ 1) ~ 1 and JZo h ( T ) dr = 1 - 
Depending on the value of Wb, the following three regimes should be distinguished: 

i) Wb <C 1. The ensemble is completely equivalent to the conventional RBM ensemble; 
diagonal matrix elements do not play an essential role; 

ii) Wi, ^> y/b. Perturbative regime. The eigenststes can be approximated by the eigen- 
states of the diagonal matrix WiSy, which are localized on single sites. The non- 
diagonal term in the Hamiltonian can be then treated via the perturbation the- 
ory; 

iii) 1 <C Wb <C \/b. Intermediate regime. It is just the regime shown to be relevant for 
the problem of two interacting particles in random potential [[|.This case is our main 
concern in the present article. 

We are going to characterize eigenfunction statistics via the following correlation func- 
tion, see 



K hm = (n\(E + 17]- H^lnYinKE -irj- 



■rn—l 

V 

l\m\ 

f, 

x exp{* ]T $t[(£ -Wt + iV A]m - i ]T //, ; <I\/.<I» ; ) (2) 

« (»i> 

where $] = (S*^ Xiii $i2i Xi 2); with S itP and Xi,p being complex commuting and Grassman- 
nian variables, respectively. The 4x4 matrices A, L are diagonal and have the following 
structure: A = diag(l, 1, —1, —1); L = diag(l, 1, —1, 1). 

Let us first calculate such a correlation function for arbitrary fixed value of the potential 
W n in the observation point n, performing both averaging over and over all Wj with 
j 7^ n (the latter averaging we denote as (...)w henceforth). Repeating all the necessary steps 
outlined in |7|] and presented in more details in || one expresses the correlation function in 
terms of the integral over the set of supermatrices Ri = T^ x P{Fi, where the supermatrices 
Pi are 4x4 block-diagonal ones and Tj belong to the graded coset space {7(1, 1/2) /?/ (1/1) x 
{7(1/ 1). The resulting expression looks as follows: 

■l—m „ 

Ki, m (E, n;r ] ) = jyjJU dRiF(R n ) exp [-i v £ Str^A - C{R)\ (3) 

i 

1 N 

C ( R ) = 9 'E(J~%StrR i R j -J2MSdet-\E - W - Ri)) w 

ij i=l 



nr{ T3 \ _ ( ^ \ ( m \ r l ~ k r m - k n k r k SdetG n 

" h \ k J V k J ( Sdet (^ -W-R, 



w 
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where G n = (E — W n — -R n ) _1 • The notations Sdet and Str stand for the graded determinant 
and graded trace, correspondingly. 

The integral over the matrices Pi can be calculated in the limit b ^> 1 by the saddle-point 
method. The saddle-point solution Pi = P s is diagonal and independent of the index i. The 
diagonal matrix elements d satisfy the following equation: 

d = /^ Wt) _J_ (4) 

Equations of similar type appeared in earlier studies of full random matrices with preferential 
diagonal [Oj, most recently in |T5,Ifi(l, and known as Pastur equation. The subsequent 



analysis depends on the value of the parameter W b characterizing the strength of the diagonal 
disorder. If Wb <C 1, we can neglect W in denominator in the r.h.s. of eq.(|^). Then the 
diagonal matrix elements distribution h(r) drops out from the formulae, and the results 
are precisely the same as for the conventional RBM ensemble. In the present paper we are 
interested in the opposite case, Wb 3> 1. Then one obviously has \d\ oc W h ~ x <C Wb and to 
the leading order in W h ~ x one finds: 



d=^P J ^ w h(W/W b )±i^h(E/W b )=Red±ilmd, (5) 



where V stands for the principal value of the integral. 

As usual, the correct saddle-point solution is equal to P s = (Red) I + i(lmd) A. In order 
to find the region of applicability of the saddle-point method we expand the functional C(R) 
around the saddle-point value and calculate the corrections due to gaussian fluctuations. 
The latter turn out to be of the order of (SPi) 2 oc b~ x . Comparing this value with the 
saddle point one P 2 ~ d 2 ~ W b ~ 2 we conclude that corrections are small as long as W 2 <C b. 
Thus, our calculation is completely legitimate everywhere in the nonperturbative regime 
1 <C Wb <C b 1 ! 2 which is just the case relevant for the physical applications of SRBM [[]]. 

Introducing the set of matrices Qi = — zT~ 1 AT i and using the identity 

! (E-W-Red)I-(lmd)Q 
(E-W n - R n ) \ P=Ps - {E _ w _ Red) 2 + [lmd) 2 

one arrives at the following expression for the correlator (^): 

K ltm (E,n',7i) = J dfi(Q)^ m (Q)e- s{Q) (6) 

where the action 

N-l N 

S(Q) = -i-Str Q^+i + ie E Str iQi A ) (7) 
^ i=i i=i 

defines the standard one-dimensional nonlinear graded a— model on a lattice characterized 
by the coupling constant 7 = (T/2) 2 J2 r J{r)r 2 and the effective level broadening e = (T/2)rj, 
the parameter T/2 being equal to T/2 = Imd cx W b ~ . 

The nonlinear a— model defined in eqs.(||[7|) was studied in much details in In 
particular, one can immediately extract the value of the localization length I which is known 
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to be proportional to the coupling constant: I = 47 oc b 2 /W 2 , in full agreement with the 
results of [§. Another quantity that can be most easily calculated is the mean local DOS 
defined in eq.(|l]) and given by: 



p(E, n) = — lm(n\(E — ir] — H) l \n) L^ 

7T 

= -lmK.i =0 , m=1 (n,E;r) -> 0) 

7T 

1 Im d 



tt(E- W n -Red) 2 + (Imd)'- 



(8) 



We conclude therefore that typically the local DOS is a Lorentzian centered around 
E = W n + Red with the width 



7T 



/ E \ 1 



r/2 = W=-^-j K - (9) 

as it was indeed found in the numerical studies [|12|]. The center of this Lorentzian is shifted 
from the local value of the random potential, W n , by the amount Red ~ Wf 1 . The shift is 
small compared to the typical values of W n ~ W^. 

Knowing all the correlators JCi m (E , n; rj) , one can extract the full set of the eigenfunc- 
tion moments P q {E) = J2 n \^a(n)\ 2q and, finally, the whole probability distribution of the 
eigenfunction amplitude |\l/ Q (n)| 2 0. Straightforwardly repeating all the necessary steps 
one finds that all moments P q {E) for SRBM are proportional to the corresponding moments 
for standard RBM at the same values of the parameters N and 7. Putting the energy E to 
be zero for the sake of simplicity, one obtains: 

p q ( E = 0) lsRBM= ([W 2 + (r/2) 2 ]«) Pq ^ E = °^ 'standard RBM ( 10 ) 

It is well known that for the standard RBM in the localized regime N 3> 7 one has P q oc r y l ~ q 
0. The relation eq.([L0|) tells us that P q IsRBM^ (T^) 1 "" K (Wg/l) 9 ' 1 . In particular, for 
the participation ratio £ = P 2 1 one has £ oc l/W 2 oc b 2 /W A , thus proving the abovemen- 
tioned difference between £ and I discovered by Jacquod and Shepelyansky ||T2"f . 

The relation between the moments eq. (0) allows one to express the distribution function 
of the normalized eigenfunction amplitude y = N\^ 2 \ for the Shepelyansky RBM model in 
terms of that for the standard RBM: 

V SRBM (V) = £ r (^-V^l) V R ™ (y T 4u) (11) 



AW b Ji (u-1) 1 / 2 \2W b ) v 4 

The actual form of the function V (y) depends on the scaling ratio iV/7 and can be 
found in Refs. |||7j It takes a simple form in both localized limit N ^> 7 and delocalized 
limit iV <C 7. For example, for the latter case V RBM (y) = e~ y and therefore one gets: 



V SRBM {y) = _ Wb _9_ 



-s/r 2 /4 roo _ 2 
dze z h 



s/y J -00 \^/yW h 



(12) 
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This distribution clearly displays the presence of two scales. All moments (y q ) = P q 
with q > 1 are dominated by the region where yT 2 ~ 1 (correspondingly, |\I/ 2 | ~ W 2 /N), 
whereas the normalization integral JV(y)dy is dominated by the values yW 2 ~ 1, where 
eigenfunction amplitude is small: |\I/ 2 | ~ (WfN) -1 . This result corresponds the following 
picture of a typical delocalized eigenstate: the eigenfunction consists of isolated peaks with 
typical amplitude l^ 2 ) ~ W b 2 /N separated by regions of a typical spatial extent L ~ W 2 
filled in with a low amplitude components |^ 2 | ~ (W 6 2 iV) _1 . 

Essentially the same picture holds for the regime of strong localization I <C N. Here 
any eigenstate has a profile which is exponentially small outside the spatial region of the 
size I. However, within this region there are isolated spikes of the amplitude |\I/ 2 | ~ W 2 /l 
separated by the low amplitude regions with a typical extent L ~ W 2 where the wave 
function amplitude is small: |\1/ 2 | ~ (H 7 " 2 /) -1 ~ 1/6 2 . 

One can also calculate for SRBM other quantities known for standard RBM. For example, 
one can be interested in level-to-level fluctuation of the participation ratio £. Performing 
such a calculation one finds that if one normalizes the inverse participation ratio by its mean 
value (P2) oc W 4 /b 2 , then the distribution of the quantity z = P/(P2) coincides exactly 
with the distribution found for standard RBM in ||. This fact suggests that envelopes of 
high amplitude peaks in SRBM are typically quite similar to envelopes of eigenfunctions in 
conventional RBM, after appropriate rescaling. 

Our last comment concerns the spectral correlator Y 2 (u) = (p(E)p(E + u)) for SRBM. 
It is easy to satisfy oneself that everywhere in the region Wb <C b 1 ^ 2 the function Yz(uS) for 
SRBM coincides with that known for RBM, as long as two-level separation u is small in 
comparison with the spectral width T of the LDOS: uj <^ W b ~ . Recently Prus and Sivan 
16| addressed the same question for the case of full matrices b = N. They found that the 



functions Y 2 (u) is given by the same expression as that for the standard Gaussian matrices 
everywhere in the region uj <ti 1/Wb as long as Wb N 1 ^ 2 . The latter result is also in 
correspondence with numerical studies by Lenz et.al [F8] who found that the crossover from 
the Wigner-Dyson statistics to the Poissonian one occurs at the scale W 2 ~ N. This scale 
is much larger than the scale Wb ~ 1 necessary to induce changes in the form of the mean 
density ||T8| , |T9f . Thus, the sparse structure of the eigenstates discussed above has no effect 
on the spectral statistics at relatively low frequency u, as long as the system stays well in 
the nonperturbative regime. 

Y.V.F is grateful to D.Shepelyansky for attracting his interest to the problem and for the 
useful communication and to O.Prus for informing him on his unpublished results This 
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